Quantum search by parallel eigenvalue adiabatic passage 
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We propose a strategy to achieve the Grover search algorithm by adiabatic passage in a very 
efficient way. An adiabatic process can be characterized by the instantaneous eigenvalues of the 
pertaining Hamiltonian, some of which form a gap. The key to the efficiency is based on the use 
of parallel eigenvalues. This allows us to obtain non-adiabatic losses which are exponentially small, 
independently of the number of items in the database in which the search is performed. 
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I. INTRODUCTION 

Quantum computation by adiabatic evolution has been 
proposed as a general method of solving search problems, 
mainly to exploit its robustness towards unitary control 
errors and decoherence [H, Q ■ In contrast to the standard 
paradigm of quantum computation which is imple- 
mented through gates embedded in a quantum circuit, 
continuous- time algorithms Q, and in particular adia- 
batic ones [H, H, @ proceed through the controlled evolu- 
tion of some Hamiltonians designed to solve the specified 
problem. The adiabatic Grover algorithm, for instance, 
involves a time-dependent Hamiltonian which smoothly 
drives the system, in a time exhibiting a quadratic 
speedup, from one of its eigenstates \ w) that is easily pre- 
pared to a connected eigenstate that coincides with the 
marked entry |m) of the database. This can be achieved 
with the two-parameter Hamiltonian 



H = a{t)Hi + b{t)H[, 



(1) 



where = \w){w\ and Hf — \m){m\ are simply projec- 
tors on the appropriate states while a and b are time- 
dependent parameters which vanish at the final and ini- 
tial times, respectively, to ensure that the prepared and 
target states are eigenstates. 

The eigenstates of H form in general an avoided cross- 
ing as a function of time. The search is achieved when the 
dynamics follows adiabatically the instantaneous eigen- 
state connected initially to the prepared state \w) and 
finally to the marked state |m). The way the parameters 
a and b vary around the avoided crossing is key to mak- 
ing the search exhibit or not a quadratic speedup. It has 
been shown that achieving a quadratic speedup requires 
a non-linear dynamics of the parameters: The dynamics 
has to slow down when approaching the smallest gap of 
the avoided crossing and has to accelerate afterwards Q . 
This strategy will be referred to as local strategy. 

In this paper we apply the strategy of optimal adi- 
abatic passage developed in Refs. [a, 0| for two-level 



models of the form (with real couplings) 



H 



A(t) fl{t) 

n{t) -A{t) 



(2) 



It has been shown that, for a given smooth pulsed-shape 
couphng of the form n{t) = iloA(t) (with A(±oo) = 0) 
and for the parametrization A(i)/Ao -I- ^l{t)/UQ = 1, the 
population transfer is the most efficient, in the adiabatic 
regime, when the instantaneous eigenvalues are parallel. 
This corresponds to Aq = Hq, i.e. to level hues (corre- 
sponding to circles of equation U{t)'^ + A{t)'^ = Ag) in 
the diagram of the difference of the eigenvalue surfaces 
as a function of the two parameters ^l{t) and A(i). 

We show that this strategy, referred to as parallel strat- 
egy, applied to the problem of quantum search using the 
two-parameter Hamiltonian ([T]) leads to a Grover type 
search, i.e. scaling with time as Vn. It is moreover 
more efficient than the local strategy proposed in Ref. 
|5] since it allows one to increase the success rate to hit 
the searched state. 

The paper is organized as follows. In Sec. II, we de- 
scribe the model and the strategies to achieve the quan- 
tum search. In particular, the local and parallel strategies 
are presented. Section HI is devoted to the definition of a 
cost and the calculation of the non-adiabatic losses which 
are used to define the optimality and to compare the lo- 
cal and parallel strategies. The comparison is illustrated 
numerically in Sec. IV while the conclusions are given in 
Sec. V. 



II. MODEL AND STRATEGIES 

The marked state \m) is one of the computational basis 
states |1), • • • , \N). The two-parameter Hamiltonian ([3]) 
can be rewritten in an orthogonal basis which features 
the marked state |m) and the uniform superposition of 
unmarked states, \u) = J^i^m TF^'*^' 
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H = 



^^-a) \n,){m\+(^ + A]\u){u\ 



+ n i\u)(m\ + \m)(u\) , 



(3) 



2 



with n = ay/N ~ 1/N and A = (a - 6)/2 - a/N. Its 
eigenvalues are 




2ab 1 



2 

N 



(4) 



and the pertaining eigenvectors read 



n 



1+) = cos0|u) +sin6'|m), tan26'=— , 6le[0,7r/2[ 

I-) =sin6'|w) -cos6l|m). (5) 

In the adiabatic representation, the Hamiltonian be- 
comes 

H^^ = B^HB-iB^B (6) 

= A+|+)(+| + A_|-)(-| + le (i-x+i - |+)(-|) , 

where the unitary transformation B is formed by the in- 
stantaneous eigenstates of i?, and with the off-diagonal 
non-adiabatic coupling 



1 - _ ^/N - 1 ab - ah 

2 A2 + f72 " ]V (A+ - A_)2 



(7) 



In the adiabatic limit, when the characteristic time of 
the process becomes arbitrarily large, the non-adiabatic 
coupling can be neglected and the dynamics follows the 
adiabatic state(s) connected to the initial state. For the 
search problem, the initial state is the uniform superpo- 
sition \w) ~ Y^^=i V) which gives no particular role 
to any state of the computational basis. Since |w) is the 
eigenvector of unit eigenvalue of iJj, taking h{t\) = im- 
plies that the instantaneous eigenvector |-|-) is connected 
to \w) at the initial time tj: 



6(ti) = 0^|+)(<i) = H. 



(8) 



At the final time tf , we require this eigenvector of higher 
eigenvalue to coincide with the marked target state, 
which is satisfied for a (if) — 0, 



a(tf) = 0-> |+)(if) = \m). 



(9) 



The adiabatic theorem [l] can be recovered from 
starting from an instantaneous eigenvector, its popula- 
tion remains larger than 1 — provided the ratio of the 
off-diagonal coupling and the gap between the eigenval- 
ues is at least e. 



A, - A_ 



(10) 



with a some multiplicative constant which fixes the en- 
ergy levels and Tunoar = tf — ti the total duration of the 
process. From (j4|) one deduces that the smallest gap is 
a/\/N while ^ yields maxt 9 = VN/Tnnear- The adi- 
abaticity condition (|10p thus implies that the computa- 
tional cost is of order N, 



aTiii, 



> 2- 



N 



(12) 



A quantum algorithm with such a linear dynamics there- 
fore does not perform better than a classical search. As 
noted in Ref. [5| , this stems from the fact that by apply- 
ing (jlOp globally, i. e., to the entire time interval, one im- 
poses a constraint on the evolution rate during the whole 
computation while the constraint is only severe where 
the gap is close to the minimum. In the next section, 
we recall the strategy proposed by Roland and Cerf Q 
which amounts to applying locally the adiabatic theorem 
for infinitesimal time intervals and adapting the rate at 
which the gap in eigenvalues is crossed. Our approach, 
which is presented afterwards in Sec. Ill C[ consists in 
following level lines on the surface of eigenvalues differ- 
ence corresponding to parallel eigenvalues, i. e., instead 
of following a given path at a varying speed, we follow 
a different path at a constant speed. This approach has 
been shown [1, 0|) in a different context, to both be ro- 
bust and strongly reduce the nonadiabatic losses. 



B. Local strategy 

The strategy proposed by Roland and Cerf Q con- 
sists in applying the adiabaticity condition locally in time 
rather than on the whole interval as in (fTO|) and corre- 
spondingly adapting the rate at which the gap A+ — A_ 
is crossed. 



A, 



A- 



(13) 



This is equivalent to fixing instantaneously the non- 
adiabatic losses to at any time applying time- 
dependent perturbation theory on the Hamiltonian 
With the parametrization a + b = a, (jH) becomes 



* 2 2 



N-1 
N 



a{l — a). 



(14) 



Hence, ([T3l) yields a differential equation for a. 



A. Linear strategy 

A naive algorithm would interpolate linearly between 
the values of a and b at the initial and final times. 



a{t) = a 
b{t) = a 



tf - t 

^linear 

t - ti 

^linear 



(11) 



ae 



N 



2 y/N^ 



N-1 , 
1 - 4 a(l - a) 



(15) 



Its implicit solution satisfying the initial condition 
a{t\) = a, which arises because of the requirement 
6(ti) — 0, reads 



a{t — ti) 



1 



1 - 2a 



A. 



A- 



(16) 



3 



At the final time, one has a{tf) = so that, denoting the 
process duration by Tiocai = if ~ ^i, one obtains 



(17) 



which shows that for a of order N''\ the search dura- 
tion scales as N^^"^, in contrast to the linear strategy for 
which, according to ([H]), it scales as N . 
The inversion of (fT6|) yields 



a{t) 



a 
2 



1 - 



sit) 



1 



with s{t) = Ti^j^- The gap A+ - A_ reads 
a 1 



A+ - A_ = 



1 



(18) 



(19) 



and its minimum is Figure [T] depicts an example of 
the dynamics of the search for a high success rate. 
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FIG. 1: (Color online) Dynamics of the local strategy search 
for = 20 and e = 1/11 (leading to oTiocai « 96). Upper 
frame: a and b (in units of 1/Tiocai); middle frame: The eigen- 
values A± (in units of l/Tiocai) exhibiting an avoided crossing; 
lower frame: Populations Pu = and 
The search is achieved with probability 0.995. 




t (units of T 



FIG. 2: (Color online) Dynamics of the parallel strategy 
search for Ft,, (t) = tanh(f/r||), iV = 20 and /Jry = 4.7. Up- 
per frame: a and b (in units of l/ry); middle frame: The 
parallel eigenvalues A± (in units of l/Tn); lower frame: Pop- 
ulations Pu = Kjijt^)!^ and Pm = |(m|<^)p. The search is 
achieved with probability 0.995. 



C. Parallel strategy 

The strategy we propose here consists in following an 
appropriate level line on the surface of eigenvalues dif- 
ference as a function of the parameters a and b of the 
Hamiltonian ([3]), corresponding to parallel eigenvalues. 
Let denote this difference where f3 is some constant 

to be chosen while, as we shall see below, the Vn arises 
to avoid energy blow up with N. From Eq. ([4]), the level 
line A+ - A_ - ^ 



is given by the ellipse 



2ab 



2 



N' 



or, in canonical form, 
4/32 



(b- 



N-1 



= 1. 



(20) 



(21) 



The initial condition ([8]), i. e. 

At the final time t[, 

^= . It follows that the parametric equation of the ellipse 



b{ti) = 0, implied that 
dSl) holds so that b{ti) = 



4 



Search cost 



bit) = h/i-^^'W 




(22) 



with F(t) a strictly monotonous function such that 
F{ti) = —1 and F{t[) = 1. Here we consider exphcitly 
the case of the (analytic) hyperbolic tangent of charac- 
teristic width Tjj, F{t) = tanh(f/T||). The eigenvalues 
(jl]) read here 



A. 




^F^t)±' 



(23) 



Figure [5] displays an example of the dynamics of 
the search, illustrating how the population transfer is 
achieved from the prepared state to the marked one. In 
the numerics, we have truncated the time interval from 
t = — 4T|| to t — 4T||. We have checked that taking a 
larger interval does not change significantly the result. 
We have chosen a situation leading to the same efficiency 
as the local strategy example shown in Fig. [TJ 

Comparing Figs. [1] and [H one also notices that the 
search is achieved in an oscillatory manner for the local 
strategy whereas it is achieved in a monotonic manner 
for the parallel strategy. This feature, as well as the en- 
hancement of the success rate for the parallel strategy, 
can be interpreted using superadiabatic basis that are 
better adapted to describe the dynamics [1, Q . In gen- 
eral, for a strategy using analytic coupling parameters (as 
considered in the parallel strategy), in the adiabatic limit, 
the non-adiabatic losses are exponentially small, i.e. of 
the form ~ g-|const.|T|| ^j^^g beyond any power of 

1/T|| while the corresponding history of the dynamics is 
smooth and monotonic. For the local strategy, the condi- 
tion (|13p prevents the analyticity of the parameters. The 
discontinuity of the coupling gives rise to losses which 
are of order 2 in the off-diagonal elements of More 
generally, a discontinuity of the n^^ derivative of the cou- 
plings corresponds to polynomial non-adiabatic losses of 
the order of (l/r"+^)^, and the corresponding dynamics 
exhibits oscillations. Indeed, the transformation leading 
to ^ can be iterated, that is one may further diagonalize 
^ , which yields higher order derivatives [§, • 

As we show more precisely in the next section, for an 
identical search cost, the parallel strategy generally en- 
hances the success rate, i. e. reduces the non-adiabatic 
losses, with respect to the local strategy, or equivalently 
reduces the cost for an identical success rate. 



In actual implementations, the time-dependent cou- 
pling parameters a{t) and b(t) can be achieved, for in- 
stance, by laser fields [ll|. As a measure of the cost 
needed to achieve the Grover search we can consider a 
quantity which is the equivalent of the total laser power. 
Note that the coupling parameters can either vary sig- 
nificantly during the whole duration of the process (e. 
g. in the local strategy) or during a small fraction only 
(e. g. in the parallel strategy). In order to account for 
both situations, we define the cost C as the product of 
the peak value flpoak of the coupling parameter a{t) and 
the effective duration of the search T^s 



C 



a 



pcak-^cfF 



(24) 



where Tcff is directly related to the characteristic time T 
of variation of a{t) and is typically a few repetitions of 
it. Toff = rT. Indeed, considering a{t) as a pulse, T is its 
characteristic width whereas T^s is the full duration for 
which a{t) is significantly different from its asymptotic 
t —^ ±oo values (hence it can be defined rigorously given 
some tolerance level) . Note that one could define the cost 
as the area under a{t) for this effective duration. The 
result would generally differ only by a numerical factor 
close to one whereas (j24p is usually more convenient to 
compute. 

For the local strategy, one deduces from ([T5|) that the 
effective duration Toff corresponds to the whole duration 



local 



tf — ti given in (|17p . The cost reads thus 



Clocal — aTiocal — 



2V]V^ 



(25) 



For the parallel strategy, we consider analytical func- 
tions such as F{t) = tanh(t/Ti|), which approach their 
asymptotic values for times t with |t| > rT^\/2. From 
(1^ . one obtains the peak value Cpoak ^ Pi^-^) 
Hence, we have 



" v/iV(7V-l) " " 



(26) 



Note that this relation seems to imply that T|| can be 
choosen arbitrarily (for instance of order N'^); However, 
the non-adiabatic losses, studied below, would then in- 
crease dramatically. 

We shall compare the two strategies at identical costs 
and then focus our attention below on the success rate of 
the search. Requiring an identical cost for both strategies 
yields 



III. COMPARISON: COST AND LOSSES 

In order to compare the local and parallel strategies, 
we shall take into account both the computational cost 
and the success rate of the search. 



2iN~lWN 
" " {N-2)e ^ ' 

For large N, the characteristic time of the squared hy- 
perbolic secant, associated with the same cost for both 



5 



strategies, is thus just {a/ fir) x Tiocai- Moreover, we can 
require the effective duration of the search to be equal for 
both strategies which simply amounts to having a = (3. 
Note that if one chooses a = /? = 1 , then the search cost 
is directly equal to the effective search duration. Since 
there is no loss of generality, we shall assume a ~ (3 = 1 
throughout. 

B. Non-adiabatic losses 

The success rate of the search is determined by the 
probability to find the system in the marked state at 



the final time. This is given by 1 — Pioss where Pioss 
corresponds to the probability of the non-adiabatic losses 
from the instantaneous eigenstate initially populated to 
the other states. These losses arise because the adiabatic 
state connected to the initial state is not strictly followed 
as 6 is not strictly zero. For the local strategy, the loss 
at the final time can be calculated exactly for any N . 
Indeed, the Hamiltonian ^ can be rewritten as 



H.. 



ad 



A+ - A- 



{|+)(+|-|->(-|+'-£(|-)(+|-|+)(-|)} 



(28) 



where A+ — A_ depends on time according to ([TO]) , and use was made of ([H]) to get A+ + A_ = a = 1. Upon 
extracting the first term which gives rise to a phase and defining a new time rit) — J 



dt one obtains a 



stationary Hamiltonian. It follows that, up to a phase, the state of the system at the new time r is thus 

mr) = 



:sin(v/l + £V)|-) + (cos(v/l + £V) ^ sin(v/l + £V) ) |+). 



(29) 



r 



From (|T4|) . we deduce that at the end of the process, 
T{tf) = iarctanViV- 1. The non-adiabatic losses are 
therefore 



-'loss, local - n SYB. 



arctan 



i.e., in the limit of large N and small e 



^loss, local ^ sin^ 



(30) 



(31) 



Note that the choice of the specific values e — l/2p with 
an integer p will give losses going to for large N . How- 
ever these choices of specific and thus non-robust values 
will not be considered here since they are outside the adi- 
abatic scope. A good measure of the losses is the upper 
boundary e^. The adiabatic regime for the local strategy 
is thus reached when e <C 1, i.e., using (fT7)) . when 



^ = - » 1 TV » 1. 



(32) 



For the level line optimization, it has been shown in 
Ref. [1] that the adiabatic regime is obtained when 
r||(A+ - A_) = 2r||/V^ > 1. Actually, for F{t) = 
tanh(</T||), we can calculate the non-adiabatic losses for 
large TV, 



loss, 11 



sech 



ttTii/? 



(33) 



This result comes from the fact that, for ^ 1, the 
model ([3]) corresponds, up to a phase, to the AUen-Ebcrly 
model ^12,]. Equation ([55]) shows the \/]V scahng of the 
search cost since taking Tu growing as ViV, 



^11 1 

= - > 1 



A> 1, 



(34) 



allows one to obtain the same arbitrarily small non- 
adiabatic losses for any A^, 



-Pi 



loss, I 







(^) 


~4exp(-^) 







(35) 



Requesting the same cost for both strategies, we de- 
duce from (P7)) . (P^ and 



er 



(36) 



The losses (|35l) . being of the form e^^^r/T^ ^re beyond any 
power of 1/7, and are thus expected to be much smaller 
than the ones given by the local strategy which are of 
order . 



IV. NUMERICAL ILLUSTRATION 

The quadratic speedup of the search using the parallel 
strategy was derived on the basis of the asymptotic result 
([55)1 obtained for large A^. We first show that the losses 
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FIG. 3: (Color online) Final losses (in logarithmic scale) as 
a function of I/7 = T\^\^/yfN (dimensionless) for A*' = 20 with 
F(t) = tanh(t/r||) (oscillating full line) its asymptotic value 
given by (|33p (non- oscillating full line, almost undistinguish- 
able from the oscillating full line), and for truncated time 
domains r = 12 and r = 8 (dashed lines). 

are also well approximated by ([33)1 for finite N . In Fig. 
[3] we plot the non-adiabatic losses obtained numerically 
as a function of I/7 defined as Ty/v^- The case N = 
20 (oscillating full line) is close to the asymptotic result 
(non-oscillating one) which, as expected from ([33]) . shows 
a strong exponential decay. The search is efficient for 
1/7 > 1> with for instance Pioss, || ~ 10~^ for 7 = 1. 
Similar results hold for other values of N . 



FIG. 4: (Color online) Final losses as a function of for e = 
1/11 for the local strategy and for the parallel strategy with 
F[i) = tanh(t/r||) with a truncated time intervals (r = 12 
and r = 8) and 7 given by Eq. (I36|l to have the same cost. 



([36l) . despite the truncation of the time interval which 
breaks the analyticity of the coupling. Note, however, 
that a truncation with too small a value for r can lead 
to a significant dependance of the losses with N . 



V. CONCLUSION 



In contrast to the local strategy, the parallel strat- 
egy uses analytical couplings on an unbounded domain. 
Hence, in practice one has to truncate this domain to 
limit the time of search. This truncation of the couplings, 
breaking the continuity, leads in general to additional 
non-adiabatic losses. In Sec. IIII Al we defined the finite 
domain through the quantity r as Tcff = rT\\ . Figure [3] 
shows the losses for r = 8 and r = 12. As expected, the 
loss becomes larger for smaller r when decreasing 7. For 
r = 12, the range of validity of the asymptotic formula 
psp is approximately I/7 ^ 3. This means that up to 
1/7^3 the additional losses due to the truncation can 
be neglected (otherwise, one could take a larger value for 
r). 



Figure H] depicts the losses as a function of N for a 
given value of e and the corresponding value of 7 given 
by (j36p with two truncations r = 8 and r ~ 12. This 
figure shows that the parallel strategy is more efficient, 
as expected from the comparison of and ([M]) with 



We have proposed a strategy to achieve the Grover 
search by adiabatic passage using parallel eigenvalues. 
We have compared this parallel strategy with the known 
local strategy which requires an adaptation of the speed 
of the dynamics with respect to the given dynamical gap 
between the eigenvalues. We have shown the superiority 
of the parallel strategy: for an identical search cost, the 
parallel strategy enhances the success rate with respect 
to the local strategy, i. e. reduces the non-adiabatic 
losses, or equivalently reduces the cost for an identical 
success rate. 

Smooth analytic coupling parameters are in princi- 
ple required for the parallel strategy. We have however 
shown numerically that a truncation of the time domain, 
which is necessary in practice, preserves the higher effi- 
ciency of the success rate of the search at identical cost. 

We have here used an hyperbolic tangent for F{t) since 
it allows one to determine analytically the population 
dynamics for large N . We have checked that other simi- 
lar shapes for F{t), for instance associated to Gaussians 
which are easily performed in the laboratory, preserves 
the advantage of the parallel strategy over the local one. 
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